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(57) Abstract 

The invention relates to a method for designing a technical system in which measurement data of a predetermined system is described 
using a substitute model. To this end, a dimension figure is determined for the quality of the substitute model by comparing the measurement 
data of the system with data determined by the substitute model. The substitute model is adapted with regard to a highest possible quality 
by using the dimension figure determined for the quality. The dimension figure is improved using additional measurement data and the 
adapted substitute model is used to design the technical system. 



(57) Zusammenfassung 

Es wird ein Verfahren zum Entwurf eines technischen Systems angegeben, bei dem MeBdaten eines vorgegebenen Systems anhand 
eines Ersatzmodells beschrieben werden. Dabei wird eine MaBzahl ftir die Gtite des Ersatzmodells ermittelt, indem die MeBdaten des 
Systems mit durch das Ersatzmodell bestimmten Daten verglichen werden. Anhand der MaBzahl fur die Gtlte wird das Ersatzmodell im 
Hinblick auf eine moglichst hohe Gtlte angepaBt. Die MaBzahl wird anhand von zusatzlichen MeBdaten verbessert und das angepaBte 
Ersatzmodell wird zum Entwurf des technischen Systems eingesetzt. 
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Beschreibung 

Verfahren und Anordnung zum Entwurf eines technischen Systems 

Die Erfindung betrifft ein Verfahren und eine Anordnung zum 
Entwurf eines technischen Systems. 

Ein Systemverhalten eines technischen Systems, z.B. einer 
verfahrenstechnischen Anlage oder eines Systems der 
GroBindustrie, hangt von zahlreichen Parametern ab. Im Rahmen 
eines Entwurf s solche eines Systems, also insbesondere beim 
Neuentwurf oder bei der Anpassung bzw. Einstellung eines 
bereits bestehenden Systems, sind Vorbedingungen, z.B. 
hinsichtlich der Wirtschaf tlichkeit oder der Umweltbelastung 
des Systems, einzuhalten. Jede Vorbedingung wird als eine 
Zielfunktion formuliert, die zu optimieren im Hinblick auf 
die anderen Zielf unktionen allgemeines Bestreben ist. 

Die Aufgabe der Erfindung besteht darin, den Entwurf eines 
technischen Systems anhand von Mefidaten eines vorgegebenen 
Systems zu ermoglichen, wobei eine gezielte Bestimmung von 
MeBdaten ein Ersatzmodell dem vorgegebenen System mit hoher 
Genauigkeit annahert. Gerade im Hinblick auf eine Optimierung 
des bestehenden Systems bzw. auf einen optimierten Neuentwurf 
eines Systems ist eine derartige Nutzung bekannter MeBdaten 
von groBer Bedeutung. 

Diese Aufgabe wird gemaB den Merkmalen der unabhangigen 
Patentanspriiche gelost. Weiterbildungen der Erfindung ergeben 
sich auch aus den abhangigen Anspruchen. 

Zur Ldsung der Aufgabe wird ein Verfahren zum Entwurf eines 
technischen Systems angegeben, bei dem Meiidaten eines 
vorgegebenen Systems anhand eines Ersatzmodells beschrieben 
werden. Es wird eine MaBzahl fur die Gute des Ersatzmodells 
ermittelt, indem die Mefidaten des vorgegebenen Systems mit 
durch das Ersatzmodell bestimmten Daten verglichen werden. 
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Anhand der Malizahl fur die Gute wird das Ersat zmodell 
dahingehend angepalit, dali es eine moglichst hohe Gute 
aufweist. Die Maiizahl fur die Gute des Ersatzmodells wird 
dadurch verbessert, dafi zusatzliche Melidaten des vorgegebenen 
5 Systems ermittelt werden. Das hinsichtlich seiner Gute 
angepalite Ersatzmodell wird zum Entwurf des technischen 
Systems eingesetzt. 

Die aus vielen unterschiedlichen realisierten Systemen 

10 vorliegenden Melidaten werden zur Beschreibung des 

Ersatzmodells eingesetzt. Mit dem Ersatzmodell wird versucht, 
das vorgegebene System moglichst gut nachzubilden. Die 
Malizahl fur die Giite der Nachbildung wird ermittelt f indem 
die realen Melidaten mit den Daten, die anhand des 

15 Ersatzmodells gewonnen werden, verglichen werden. Eine grolie 
Differenz zwischen den Melidaten und den Daten des 
Ersatzmodells entspricht einer schlechten Giite, also einer 
schlechten Abbildung des vorgegebenen Systems in das 
Ersatzmodell. Durch die Malizahl fur die Gttte wird das 

20 Ersatzmodell dahingehend angepafit, dali die Gute selbst 
moglichst hoch wird und somit das Ersatzmodell das 
vorgegebene System moglichst gut beschreibt. Gibt es einzelne 
Melidaten, die von den Daten des Ersatzmodells besonders 
abweichen, so werden zusatzliche Melidaten in einer 

25 vorgegebenen Umgebung dieser Melidaten ermittelt, so dafi eine 
genauere Anpassung (und damit eine weitere Erhohung der Giite) 
des Ersatzmodells an das vorgegebene System erfolgen kann. 
Das so gewonnene Ersatzmodell mit hoher Giite wird zum Entwurf 
des technischen Systems eingesetzt. 

30 

Unter Entwurf werden allgemein verstanden sowohl der 
Neuentwurf eines technischen Systems als auch die Anpassung 
bzw. Adaption eines bereits vorhandenen technischen Systems. 

35 Eine Weiterbildung besteht darin, dali das Ersatzmodell ein 
Regressionsmodell ist. Das Regressionsmodell geht von der 
Beschreibung 
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Yi = ^j(*i) + ei 



aus, wobei 
5 (yi, Xi) vorgegebene Wertepaare, 

fp eine Funktion, die von einem Parameter p abhangt 

und 

ei einen Fehler 

bezeichnen. 



Nun ist der Fehler (als Funktion von p) zu minimieren : 



Se? = <p(p). 
i = l 

15 Geht man von folgendem Beispiel 

y = Po + Pi* + P2* 2 + e 

aus, so ist der funktionale Zusammenhang quadratischer 
20 Ordnung, das Regressionsmodell (Funktion, abhangig von P) 
hingegen ist linear. 

Eine zusatzliche Erhohung der Gute des Ersatzmodells wird 
dadurch erreicht, dafi mindestens ein Teilbereich des 
25 Ersatzmodells lokal angepaflt wird. Hierbei erfolgt eine 
Anpassung des Ersatzmodells durch zusatzliche (lokale) 
Ermittlung von Mefiwerten, die zur genaueren Anpassung des 
Ersatzmodells an das vorgegebene System herangezogen werden. 

30 Die Gute kann in einer anderen Weiterbildung anhand einer 
quadratischen Abweichung der MeBdaten von den durch das 
Ersatzmodell bestimmten Daten ermittelt werden. Die Anpassung 
des Ersatzmodells erfolgt durch Minimierung der quadratischen 
Abweichung. 

35 
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Eine Ausgestaltung besteht darin, daft die Mefidaten nach ihrer 
Giite, bezogen auf deren Abweichung von den durch das 
Ersatzmodell bestimmten Daten, sortiert werden und eine 
vorgegebene Anzahl von n% schlechtesten Mefidaten aussortiert 
5 werden. Es wird also eine Giite fur jedes Mefidatum bestimmt, 
wobei die Menge der Mefidaten, vorzugsweise in Form einer 
Liste, nach ihrer Giite sortiert und die n% schlechtesten bzw. 
die n schlechtesten Mefidaten aussortiert werden. Insbesondere 
ist zu priifen, ob die n% bzw. die n schlechtesten Mefidaten in 
10 einem zusammenhangenden Bereich liegen. Ist dies der Fall, so 
werden diese Mefidaten nicht aussortiert, da sie mit hoher 
Wahrscheinlichkeit keine Mefifehler, sondern einen 
zusammenhangenden Bereich bestimmen, der durch das 
Ersatzmodell nicht ausreichend genau abgebildet wurde . 

15 

Eine andere Weiterbildung besteht darin, dafi die Mefidaten 
einer Vorverarbeitung unterzogen werden. Da in einem realen 
vorgegebenes System eine grofie Menge an Mefidaten pro 
Zeiteinheit anf alien, ist es sinnvoll, diese Mefidaten einer 
20 Vorverarbeitung zu unterziehen und somit zu gewahrleisten, 
dafi weitgehend signifikante Mefidaten in die Bildung des 
Ersatzmodells einfliefien. Bevorzugt findet die 
Vorverarbeitung in einer Reduktion der Anzahl von Mefidaten 
ihre Auspragung. 

25 

Dabei werden die Mefidaten entsprechend vorgegebener Kriterien 
in Klassen eingeteilt. Die Mefiwerte einer Klasse werden 
bewertet und diejenigen Mefiwerte, deren Bewertung unterhalb 
eines vorgegebenen ersten Schwellwertes liegt, werden 

30 aussortiert. Durch die Aussortierung der Mefiwerte ergibt sich 
eine Reduktion hinsichtlich der Anzahl der Mefiwerte. Somit 
liegen fiir eine Weiterverarbeitung eine deutlich reduzierte 
Anzahl von Mefiwerten vor. Die Weiterverarbeitung kann mit 
gegenuber der nicht reduzierten Anzahl von Mefiwerten 

35 geringerem Rechenauf wand erfolgen. 
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Auch konnen die Klassen selbst bewertet werden. Insbesondere 
kann eine Klasse, deren Bewertung unterhalb eines 
vorgegebenen zweiten Schwellwerts liegt, aussortiert werden. 
Hierdurch ergibt sich eine zusatzliche Reduktion der Anzahl 
der Meliwerte. 

Eine andere Weiterbildung der Vorverarbeitung besteht darin, 
daft ein Kriterium fur die Klasseneinteilung darin besteht, 
daft pro Klasse Meliwerte zu einer Vorgabe von 
Einstellparametern des technischen Systems bestimmt werden. 
Typischerweise wird das technische System anhand einer 
vorgegebenen Anzahl von Einstellparametern eingestellt, nach 
Einstellung erfolgt eine (zumeist zeitverzogerte) Reaktion 
des Systems auf die Einstellparameter (Einschwingverhalten, 
Einschwingvorgang des Systems) Nach Einstellung. werden somit 
eine bestimmte Menge an Mefiwerten aufgenommen, die dem 
Einschwingvorgang zugeordnet werden konnen, wobei nach 
abgeschlossenem Einschwingvorgang (Obergang in den 
stationaren Betrieb) weiterhin Meliwerte anfallen, die dem 
vorgegebenen Satz Einstellparameter zugeordnet werden. Durch 
Verstellung der Einstellparameter wird eine neue Klasse 
bestimmt. Alle Meliwerte, die jeweils nach einer Verstellung 
der Einstellparameter anfallen, gehoren in eine eigene 
Klasse. 

Zusatzlich kSnnen Meliwerte einer Klasse, die dem jeweiligen 
Einschwingvorgang zuordenbar sind, aussortiert werden. 
Weiterhin konnen fehlerhafte Meliwerte aussortiert werden. Die 
Einstellung groJier technischer Systeme ist in vielen Fallen 
auf eine langf ristigen stationaren Betrieb ausgerichtet . 
Meliwerte, die sich auf den Einschwingvorgang (von kurzer 
Dauer im Verhaltnis zum stationaren Betrieb nach 
abgeschlossenem Einschwingvorgang) beziehen, werden sinnvoll 
aussortiert, da durch sie Meliwerte fur den stationaren 
Betrieb verfalscht werden. Insbesondere im Rahmen einer 
Modellierung des technischen Systems, sind die Meiidaten des 
stationaren Verhaltens des technischen Systems interessant. 
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Eine Ausgestaltung besteht darin, die Anzahl der Mefiwerte in 
einer Klasse dadurch zu reduzieren, dafi mindestens ein 
reprasentativer Wert fur die Mefiwerte der Klasse bestimmt 
wird. Solch ein reprasentativer Wert kann sein: 

a) ein Mittelwert (z.B. ein gleitender Mittelwert) der 
Mefiwerte der Klasse, 

b) ein Maximalwert der Mefiwerte der Klasse, 

c) ein Minimalwert der Mefiwerte der Klasse, 

d) ein Median. 

Bei Variante d) liegt ein Vorteil darin, dafi immer ein Wert 
bestimmbar ist, den es tatsachlich gibt, wohingegen der 
Mittelwert a) selbst nicht als Wert vorkommt. 

Je nach Anwendungsf all, kann eine geeignete Wahl zur 
Bestimmung des reprasentativen Werts einer Klasse erfolgen. 

Eine ganze Klasse mit Mefiwerten kann aussortiert werden, wenn 
diese weniger als eine vorgegebene Anzahl Mefiwerte enthalt. 

Eine andere Ausgestaltung besteht darin, dafi diejenigen 
Mefiwerte aussortiert werden, die urn mehr als einen 
vorgegebenen Schwellwert von einem vorgebbaren Wert 
verschieden sind. Der vorgebbare Wert kann ein Mittelwert 
aller Mefiwerte der Klasse oder ein zu erwartender Mefiwert auf 
die jeweiligen Einstellparameter des technischen Systems 
sein. 

Im Rahmen einer anderen Weiterbildung werden die mittels 
Entwurf gewonnenen Daten zur Steuerung einer technischen 
Anlage eingesetzt. Zusatzlich kann die Steuerung der 
technischen Anlage zur Laufzeit des Systems, also Online, 
erf olgen. 

Auch wird zur Losung der Aufgabe eine Anordnung zum Entwurf 
eines technischen Systems angegeben, die eine 
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Prozessoreinheit aufweist, welche Prozessoreinheit derart 
eingerichtet ist, dafi Mefidaten eines vorgegebenen Systems 
anhand eines Ersatzmodells beschreibbar sind. Eine Mafizahl 
fur die Gute des Ersatzmodells ist ermittelbar, indem die 
5 Mefidaten des vorgegebenen Systems mit durch das Ersatzmodell 
bestimmten Daten verglichen werden. Aus der Mafizahl fur die 
Gute ist das Ersatzmodell dahingehend anpafibar, daft es eine 
moglichst hohe Gtite aufweist. Die Mafizahl fiir die Gute des 
Ersatzmodells ist weiter verbesserbar, indem zusatzlich 
10 Mefidaten des vorgegebenen Systems ermittelt werden. Das 
hinsichtlich seiner Gute angepaflte Ersatzmodell ist zum 
Entwurf des technischen Systems einsetzbar. 

Diese Anordnung ist insbesondere geeignet zur Durchfuhrung 
15 des erf indungsgemafien Verfahrens oder einer seiner vorstehend 
erlauterten Weiterbildungen. 

Ausfuhrungsbeispiele der Erfindung werden nachfolgend anhand 
der Zeichnung dargestellt und erlautert. 

20 

Es zeigen 

Fig.l ein Blockdiagramm, das Schritte eines Verfahrens zum 
Entwurf eines technischen Systems enthalt; 

25 

Fig. 2 eine schematische Skizze eines Recovery-Boilers; 

Fig. 3-5 Eingangsgrofien, Stellgrofien und Ausgangsgrofien des 
Recovery-Boilers . 

30 

In Fig.l ist ein Blockdiagramm dargestellt, das Schritte 
eines Verfahrens zum Entwurf eines technischen Systems 
enthalt. In einem Schritt 101 wird anhand von Mefidaten ein 
35 Ersatzmodell gebildet. Bevorzugt ist dieses Ersatzmodell ein 
Regressionsmodell . Dann sind bei einer quadratischen Funktion 
f die erste und die zweite Ableitung unmittelbar gegeben. Urn 
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clas in Schritt 101 entstandene Ersat zmodell auf die Mefidaten 
anzupassen, also eine Verfeinerung des Ersatzmodells 
vorzunehmen, so dafi die Mefidaten das Ersatzmodell in 
ausreichender N&herung beschreiben, wird in einem Schritt 102 
5 eine Mafizahl fur die Giite des Ersatzmodells bestimmt. Diese 
Mafizahl wird bestimmt, indent die Mefidaten des vorgegebenen 
Systems mit durch das Ersatzmodell bestimmten Daten 
verglichen werden. Bevorzugt erhalt jedes Mefldatum eine 
Mafizahl fiir die Giite, welche Mafizahl die Abweichung des 

10 Mefidatums mit dem zugehorigen von dem Ersatzmodell bestimmten 
Wert kennzeichnet . Die Summe aller Mafizahlen fiir die Giite fur 
alle Mefidaten bestimmt eine Gesamt giite fiir das Ersatzmodell. 
In einem Schritt 103 wird die Giite maximiert, indem die 
Mafizahl fur die Giite maximiert wird. Zusatzlich wird das 

15 Ersatzmodell verfeinert, also die Giite der Obereinstimmung 
des Ersatzmodells mit dem vorgegebenen System erhoht, indem 
gezielt Mefidaten an den Stellen bestimmt werden und in die 
Bildung des Ersatzmodells einfliefien, an denen eine grofie 
Abweichung zwischen Mefidaten und durch das Ersatzmodell 

20 bestimmten Daten besteht. 1st eine entsprechend hohe Giite fiir 
das Ersatzmodell bestimmt, so wird in einem Schritt 104 
dieses Ersatzmodell fiir den Entwurf des technischen Systems 
eingesetzt. Der Entwurf kann sowohl ein Neuentwurf (vgl. 
Schritt 105) oder eine Anpassung eines schon bestehenden 

25 technischen Systems (vgl. Schritt 106) sein. 

Fig. 2 zeigt eine schematische Skizze eines Recovery-Boilers. 
Nachfolgend wird anhand des Beispiels "Recovery-Boiler" ein 
Ausfiihrungsbeispiel des oben beschriebenen Verfahrens 
30 veranschaulicht . 

In der Papier- und Zellstof f industrie werden zum Aufschlufi 
von Zellstoff verschiedene Chemikalien sowie Warme und 
Elektroenergie benotigt. Aus einer eingedeckten Prozefiablauge 
35 (Schwarzlauge) lassen sich mit Hilfe des Recovery-Boilers die 
verwendeten Chemikalien und zusatzlich Warmeenergie 
zuruckgewinnen. Ein Grad fur die Zuruckgewinnung der 
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Chemikalien ist von entscheidender Bedeutung fur die 
Wirtschaf tlichkeit der Gesamtanlage. 

Die Schwarzlauge wird in einem Schmelzbett 201 verbrannt. 
5 Dabei bildet sich eine Alkalischmelze, die iiber eine Leitung 
202 abfliefit. Aus den Bestandteilen der Alkalischmelze werden 
in weiteren Verf ahrensschritten die eingesetzten Chemikalien 
zuruckgewonnen. Frei gewordene Verbrennungswarme wird zur 
Erzeugung von Wasserdampf genutzt. Die Verbrennung der 

10 Ablauge und damit die Ruckgewinnung der Chemikalien beginnt 

mit der Zerstaubung der Schwarzlauge iiber Zerstauberdusen 204 
in eine Brennkammer 203. Partikel der zerstaubten 
Schwarzlauge werden bei ihrem Fall durch das heifie Rauchgas 
getrocknet. Die getrockneten Laugenpartikel fallen auf das 

15 Schmelzbett 201 , wobei eine erste Verbrennung und eine 

chemische Reduktion stattfinden. Fluchtige Bestandteile und 
Reaktionsprodukte gelangen in eine Oxidationszone, in der 
oxidierende Reaktionen ablaufen und in der die Verbrennung 
abgeschlossen wird. 

20 

Wichtige Zielvorgaben fur die Steuerung des Recovery-Boilers 
sind die Dampf produktion zur Energiegewinnung, die Einhaltung 
von Emissionswerten unter Umweltgesichtspunkten und die 
Effizienz der chemischen Reduktion. 



Der Verbrennungsvorgang, und damit die Zielvorgaben, werden 
insbesondere durch die Luftzufuhr in drei Ebenen (Primary Air 
(PA), Secondary Air (SA) , Tertiary Air (TA) ) gesteuert. Der 
Gesamtprozefi unterliegt zahlreichen Einflussen, die bei der 
30 Modellierung zu berucksichtigen sind: 



25 



a) Die Messung der Grolien unterliegen oftmals starken 
Schwankungen ; 



35 



b) Es existieren nicht gemessene und nicht mefibare 
Einf luligroiien; 
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c) Bei jecier Anderung der einstellbaren Parameter komnit 
es zu Einschwingvorgangen; 

d) Die technische Anlage verschmutzt und wird in 
vorgegebenen Abstanden gereinigt, wodurch im 
Systemverhalten jeweils ein zeitlicher Drift bewirkt 
wird. 

Die gemessenen Groflen des Gesamtprozesses werden in 
Eingangsgroflen (vgl. Fig. 3) und Ausgangsgroflen (vgl. Fig. 5) 
unterteilt. Jede Minute werden Meflwerte abgespeichert . Vier 
der Eingangsgroflen sind gleichzeitig auch Stellgroflen (auch: 
einstellbare Parameter; vgl. Fig. 4). Die Stellgroflen sind im 
wesentlichen als unabhangig voneinander einstellbare freie 
Parameter des Gesamtprozesses anzusehen. Einige der anderen 
Eingangsgroflen sind von den Stellgroflen mehr oder minder 
abhangig. Gemafl einer Vorgabe sind beim Recovery-Boiler die 
Groflen "BL Front Pressure" und "BL Back Pressure" stets 
gleich zu regeln. Die vier Stellgroflen (vgl. Fig. 4) sind 
vorzugsweise abzuspeichern als Stellgroflen (mit dem 
gewlinschten, voreingestellten Wert) und als Eingangsgroflen 
(mit dem gemessenen, realen Wert) . 

Beim Recovery-Boiler besteht eine Problemstellung darin, in 
Abhangigkeit von den einstellbaren Parametern bestimmte 
Zielvorgaben, die uber gemessene Groflen definiert werden, zu 
erfullen. Hier wird eine dreistufige Vorgehensweise zur 
Losung des Problems gewahlt: 

1. Die zu betrachtenden Zielvorgaben werden durch 

stochastische Methoden modelliert, wobei diese Modelle 
durch neue Messungen aktualisiert werden 

(datengetriebene, empirische Modellierung) . Dabei ist es 
sinnvoll, nicht nur ein einziges Modell zu verwenden, 
sondern globale Modelle fur die Identif ikation 
interessanter Gebiete in einem durch die Zielvorgaben 
bestimmten Parameterraum und lokale Modelle zur exakten 
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Berechnung optimaler Arbeitspunkte einzusetzen. Die 
verwendeten Modelle warden durch Giitemalie bewertet. 

2. Falls die betrachteten Modelle aufgrund der Datenlage 
5 nicht hinreichend genau sind (Gutemafle) , werden gezielt 



3. Fur die lokale Optimierung werden lokale Modelle 

konstruiert und die zur Verfiigung stehenden Datensatze 
gegebenenfalls gezielt erweitert (Experimental Design) . 

15 

Bei den Zielvorgaben handelt es sich um physikalisch- 
technische bzw. betriebswirtschaf tliche Kriterien, die in der 
Regel Randbedingungen und/oder Sicherheitsbedingungen 
entsprechen mussen. Haufig sind mehrere dieser Kriterien 

20 gleichzeitig zu betrachten. Die Verwendung eines 

stochastischen Modells kann insbesondere dazu verwendet 
werden, die zu optimierenden ZielgroBen und ihre Abhangigkeit 
von den einzustellenden Parametern irn Rechner zu simulieren. 
Dies ist dann notwendig, wenn Messungen sehr kostenintensiv 

25 bzw. sehr zeitauf wendig sind. Bei Sicherheitsanf orderungen 
konnen mogliche Gef ahrensituationen vermieden werden. 

Beim Recovery-Boiler ist eine Online-Optimierung, die auf 
mehreren Daten basiert, notwendig, weil die physikalisch- 

30 chemischen Prozesse nicht mit ausreichender Genauigkeit 

quantitativ modelliert werden konnen und weil das Verhalten 
der Anlage im Verlauf des Betriebs Schwankungen unterliegt. 
Das Wissen uber dieses Verhalten muli stetig durch gezielte 
Wahl neuer Arbeitspunkte erweitert werden. Daher empfiehlt 

35 sich im Rahmen der Online-Optimierung das bereits 

beschriebene dreistufige Vorgehen der stochastischen 
Modellierung und der mathematischen Optimierung. 



10 



neue Arbeitspunkte zur Modellverbesserung ausgewertet 
(Experimental Design) . Ferner werden durch Verwendung 
globaler stochastischer Optimierungsverf ahren bzgl. der 
Zielvorgaben attraktive Gebiete in Abhangigkeit vom 
aktuellen globalen Modell identif iziert . 
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BESCHREIBUNG DER EINGANGSGROSSEN 

5 Die a EingangsgrSfien (a € N, N: Menge der naturlichen 

Zahlen) sind im allgemeinen von n Stellgroflen n e N und von 
Zuf allsef fekten abhangig. Sie konnen wie folgt beschrieben 
werden: 

10 Es seien (£l,S,P) ein Wahrscheinlichkeitsraum und S 7 eine 
Borelsche a-Algebra uber R V (R: Menge der reellen Zahlen) 
ftir jedes v e N . Die Eingangsgrofien werden uber eine 
^ n x S - 2* a -mefJbare Abbildung cp dargestellt: 

15 cp : R n x Q -> R a (1) m 

Die Def initionsmenge der Abbildung <p ist ein kartesisches 
Produkt zweier Mengen. Betrachtet man die jeweiligen 
Projektionen auf die Einzelmengen, so erhalt man folgende 
20 Abbildungen : 

<p x :Q -> R a , co -» <p{x,(o) fiir alle x e R n (2), 
tp^iR 11 R a , x -» <p(x, a>) fiir alle co e Q (3) . 

25 

{<p x ; x e R n | ist ein stochastischer Prozefi mit einer 

Indexmenge R n und eine Abbildung cp® ist fiir jedes Ereignis 
© e Q ein Pfad dieses stochastischen Prozesses. 

30 Beim Recovery-Boiler ist n=4 und a~14 (nach Elimination der 
Gr6fle "BL Back Pressure"). 

Aufgrund der geforderten Meftbarkeit der Abbildung cp x ist fur 
jedes x € R n die Abbildung q> x eine Zuf allsvariable . Unter 
35 geeigneten Zusatzvoraussetzungen konnen Erwartungswerte und 
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hohere Momente betrachtet werden. Dieser Zugang ermoglicht 
den Schritt von stochastischen Modellen zu deterministischen 
Optimierungsproblemen . Bei einem deterministischen 
Optimierungsproblem ist die Zielfunktion direkt mittels einer 
5 Variablen einstellbar, wohingegen die stochastische Grolie die 
Zielfunktion beeinflufit, aber keine gezielte Einstellung 
ermoglicht . 

10 BESCHREIBUNG DER AUSGANGSGROfiEN 

Das Prozeflmodell M des Recovery-Boilers wird als Funktion in 
Abhangigkeit von den Eingangsgroften und weiteren 
Zuf allsef f ekten beschrieben. Dabei sei (Q,S, P) der obige 
15 Wahrscheinlichkeitsraum. Das Prozeftmodell M ist dann eine 
& a x S - ? b -me(5bare Abbildung: 

M : R a x Q R b ( 4 ) , 

20 wobei b die Anzahl der Ausgangsgroften bezeichnet. 

Da der Recovery-Boiler einem zyklischen zeitlichen Drift 
unterliegt (von Reinigungsphase zu Reinigungsphase) , ist 
zudem eine Beschreibung mit einem Zeitparameter denkbar. Die 
25 Ausgangsgrolien lassen sich durch 2* n x s - ^ b -meftbare 
Abbildungen v|/ darstellen: 

y : R n x CI -» R b (5) , 

30 (x, to) m(<p(x, co), co) ( 6) . 

Betrachtet man die jeweiligen Projektionen auf die 
Einzelmengen der Def initionsmenge, so erhalt man folgende 
Abbildungen 

35 

V|/ x : Q R b , Q) i|/(x, <o) flir alle x e R n (7), 
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xj/^R" -> R b , x -» \|/(x, <d) fiir alle o> 6 £2 (8). 

{v x ; x e R n } ist ein stochastischer Prozeft mit einer 

Indexmenge R n und die Abbildung vj/® ist fiir jedes co € Q ein 
Pfad dieses stochastischen Prozesses. 

Beim Recovery-Boiler ist b=15. 

Die Tatsache, dafi bei der Definition von vj/ zwischen den 
verwendeten Ereignissen cd nicht unterschieden wird, bedeutet 
keine Einschrankung, da CI als kartesisches Produkt aus einem 
Ql und einem CI2 dargestellt werden kann. Die obige 
Darstellung umfaBt somit auch das Modell: 

\|/ : R n x Cl ± x Q 2 -> R b 
(x, ©i, ©2) -» M(q>(x, coi), CD2) 



(9) , 

(10) . 



BESCHREIBUNG DER ZUR VERFOGUNG STEHENDEN DATENSATZE 

Mit den Beschreibungen in den beiden vorangegangenen 
Abschnitten kann man die Eingangsgrofien und die Ausgangs- 
groiJen gemeinsam zu MeJJgrofien (=MefJdaten) <D zusammenf assen . 
& ist eine 2> n x S - ff m -mefibare Abbildung mit m = a + b und 

O : R n x Q -» R m 

Betrachtet man wieder die jeweiligen Projektionen auf die 
Einzelmengen der Def initionsmenge, so erhalt man folgende 
Abbildungen : 
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O x : n -»• R m , a> <D(x, a>) ftir alle x e R n (13), 
5 tD^tR" R m , x -> <D(x, <o) fxlr alle co e Q (14) . 

{<D X ; x e R n } ist ein stochastischer Prozefl mit einer 

Indexmenge -R n und die Abbildung ist fur jedes to e Q ein 

Pfad dieses stochastischen Prozesses. 

10 

Fur jedes gewahlte Stellgrofientupel x werden beim Recovery- 
Boiler viele Realisierungen von O x ermittelt und 

abgespeichert, d.h. zu jedem Xj e R n werden zahlreiche 

Realisierungen 

15 

<D jk : = 0>(xj, ©j k ) (15) 

mit ©jk e Q; k = 1,2, . . . , v j ; 

20 v j e N ; j = 1,2, . . . , u ; u e N 

betrachtet. Die gespeicherten Datensatze Dj k des Recovery- 
Boilers sind also (n + m)-Tupel: 

25 D j k = (<Dj j J ' k = 1,2, . . . , vj ; j = 1,2, . . . , u ( 16) . 

Dabei wird Dj lkl vor Dj 2 k 2 abgespeichert, wenn 
Ol < 32) v ((ji = j 2 ) a (k! < k 2 )) 

30 

gilt. 



DAT EN KOM P RE S S I ON DURCH KLASSENEINTEILUNG DER PARAMETER 
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Da fur jedes Stellgrdfientupel x im allgemeinen mehrere 
Realisierungen von <I> X vorliegen, bietet sich aufgrund der 
komplexen stochastischen Eigenschaf ten des zu betrachtenden 
Prozesses als erster Schritt der statistischen Datenanalyse 
eine Klasseneinteilung der Parameter durch Bildung 
arithmetischer Mittelwerte an. Zudem werden of f ensichtlich 
fehlerhafte Datensatze ausgesondert . Ein of f ensichtlich 
fehlerhafter Datensatz ist bspw. eine physikalisch unmogliche 
Messung, der insbesondere aufgrund einer vorgenommenen 
Einstellung real gar nicht vorkommen kann. 

Vorgehensweise : 

1. Datensatze, bei denen die GrofJe "BL Front Pressure" 
ungleich der Grofie "BL Back Pressure" ist, werden 
aussortiert, da diese beiden Werte nach Vorgabe der 
Anlagensteuerung gleich sein miissen. Der Datenverlust ist 
sehr gering. 

2. Die Datensatze werden auf Klassen aufgeteilt, in denen 

die vier Einstellparameter (PA, SA, TA, BL Front 

Pressure, siehe oben) zeitlich auf einanderf olgend 

konstant sind, d.h. die j-te Klasse besteht aus den 
Datensatzen Dj # . 

3. Klassen, in denen sich weniger als 30 Datensatze 
befinden, werden ausgesondert, damit Einschwingvorgange 
keinen groBen Einfluft haben. 

4. Fur jede Klasse werden ein arithmetischer Mittelwert 
und eine empirische Standardabweichung sj fur alle 
Mefigrofien ermittelt: 
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Vj -1 



v j 



1^ 



k=l 



k = l 



(18) 



10 



5. Klassen bei denen die Mittelwerte fur die Grofien PA f SA, 
TA Oder BL Front Pressure zu weit von den entsprechenden 
Einstellparametern entfernt sind, werden ausgesondert . 
Auf diesen Klassen konnten also die Einstellwerte nicht 
erreicht werden. 



STATISTISCHE KENNGROftEN FOR DIE GEGEBENEN KLASSEN UND IHRE 
GRAPHISCHE DARSTELLUNG 

Neben den arithmetischen Mittelwerten und den empirischen 
15 Standardabweichungen, die fur die einzelnen Klassen bestimmt 
wurden, wird noch eine gemeinsame Standardabweichung s 
bestimmt gemafi 



20 



s = 



j = l 



^2 



D = l 



(m) 2 
j 



1 



(19) 



Dabei steht u fur die Anzahl der Klassen (hier 205) und v fur 
die Summe der v j , d.h. v ist die Anzahl aller verwendeten 

Mefiwerte (hier 38915) . 
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LINE ARE REGRESS IONSMODELLE FUR FUNKTIONSAPPROXIMATIONEN 

Fur jede Mefigrofce (Mefidatum) (i=l, 2, . . . ,m) wird, 

5 basierend auf dem arithmetischen Mittel uber die Klassen, ein 
lineares Regressionsmodell in Abhangigkeit von der 
quadratischen Kombination der vier Einstellparameter 
berechnet. In der folgenden Darstellung ist x e R 4 , wobei 

10 x^ : Primary Air (PA) 

(2) 

x v 7 : Secondary Air (SA) 



15 



(3) 

x v ' : Tertxary Air (TA) 



x^ : Black Liquor (BL) Front Pressure 



gilt- u € N bezeichnet die Anzahl der Klassen. Jede Mefigrofie 
20 wird durch 

<D^(x, co) = ajr(x) + e±(&) (20) 



25 



30 



mit aj_ e R 15 modelliert . Dabei gelten 

r : R 4 R 15 (21) 

fa, ft, ft, ftf -+ ( If ft, ft, 53^ ft' Sir & S|' ^ 

ftft, ftft , ftft, ftft, ftft, ftft ) T ( 22 > , 

d.h. Polynome zweiten Grades werden an die Meftdaten angepafit, 
und 

e ± : Q -> R (23) 



35 
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ist eine Zuf allsvariable mit Erwartungswert 0. 

Der Vektor a± wird mit der Methode der kleinsten 
Fehlerquadrate bestimmt, allerdings werden statt der 

Originaldatensatze (xj, die arithmetischen Mittel 

( x j' ^D^) verwendet. Diese Vorgehensweise eignet sich, da 

durch lineare Regressionsmodelle insbesondere Erwartungswerte 
geschatzt werden. Somit erhalt man das Minimierungsproblem: 



mm 



as eR 



15 



r(xif 
r ( X u) T ; 



2 J 



(24) 



Es sei ai der optimale Punkt des quadratischen 
Minimierungsproblems aus Gleichung (24). Ferner gelte: 









( # 1 








-(15) 

^ a i J 



6 R 



u 



(25) , 



u 



yi : = 



U 3 = 1 



e R 



(26) 



Zur Validierung des Regressionsansatzes wird ein 
Bestimmtheitsmafi R 2 berechnet gemafi 



FT : = 



yjyi - uyf 



^(i)T^(i) _ uy 2 



(27) 



mit 
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(28) 



Je naher R£ bei 1 liegt, desto besser wird die abhangige 
Variable durch die unabhangigen Variablen dargestellt 
(o < R? < l). 

Zusatzlich wird ein Maximum E^x fur einen Absolutwertes der 
Abweichung der Daten vom Modell angegeben als 



10 



,(i) 
j max 



max 
j = l,...,u 



(29) 
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E SK)% ist derjenige Wert, unterhalb dessen mindestens 90% der 

Absolutwerte der Abweichungen der Daten vom Modell liegen. 

Analog dazu ist Eg*Q% derjenige Wert, unterhalb dessen 

mindestens 80% der Absolutwerte der Abweichungen der Daten 
vom Modell liegen. Mit dem optimalen Punkt a^ des 

Minimierungsproblems gemali Gleichung (24) lafit sich ein 

Modell des Erwartungswertes der Meftgrofie angeben zu 



20 



= R 



n 



R 



(30) , 



afr(x) 



(31) 



Insbesondere lalit sich der Gradient VO^ analytisch angeben 
25 mit 



VdCD(x) = ^£ ( x ) . a fur alle x e R n 
dx 



(32) 



30 



OPTIMIERUNGSKONZEPTE UND DATENGETRIEBENE MODELLANPASSUNG 
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ALLGEMEINE VORGEHENSWEISE 



Die Steuerung des Recovery-Boilers erfolgt durch einen 
5 Anlagen-Operator . Aufgrund von Vorgaben des Anlagen- 
Operators, ergibt sich eine Funktion 



F : R n x fi -» R (33) , 

10 (x, co) -> g(o(x, co), x) (34). 

In die Wahl der Funktion 



15 



g : R m x R n -+ R 
(r\, x) -» g(ti, x) 



(35) , 
(36) 



gehen die Vorgaben des Anlagen-Operators ein. Ein Ziel 
besteht in der globalen Minimierung des Erwartungswertes der 
20 Funktion F. Aufgrund der stochastischen Komplexitat des 
beschriebenen Problems bietet sich eine dreistufige 
approximative Vorgehensweise an: 

Mit Hilfe der oben dargestellten stochastischen 
Beschreibung der Mefidaten durchgef uhrten statistischen 
Datenanalyse und der daraus resultierenden Modelle O fiir 
den Erwartungswert von <D , wird das folgende globale 
Modell F fiir den Erwartungswert der Funktion F bestimmt: 



1. 

25 



30 



F : R n — ► R 



x -» g(c»(x), x) 



(37) , 

(38) . 



35 



Fur konvexe Funktionen g kann man mit Hilfe der 
Jensenschen Ungleichung die folgende Beziehung angeben: 
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F(x) <; e(f(x,#)) (39) , 

d.h. F approximiert den Erwartungswert von F von unten. 
Das Symbol "•" deutet an, daft die Zuf allsvariable durch 
die Beziehung E(...) eliminiert wird. 

2. Unter Verwendung des entwickelten globalen Modells konnen 
Punkte in geeigneten Umgebungen globaler Minimierer von 

F durch stochastische Methoden zur globalen Minimierung 
(siehe [1]) berechnet werden. Diese Punkte dienen als 
Startpunkte fur eine lokale Minimierung. Falls das 
betrachtete globale Modell aufgrund der Datenlage nicht 
hinreichend genau ist, also die Gute nicht ausreichen 
hoch ist, werden gezielt neue Arbeitspunkte zur 
Modellverbesserung ausgewertet (Experimental Design) . 

3. Fur die lokale Minimierung ist das globale Modell F 
eventuell nicht genau genug. Daher wird jeweils ein 
lokales Modell des Erwartungswertes von F verwendet. Ftir 
die Modellierung verwendet man wieder geeignete lineare 
Regressionsmodelle, wobei jetzt bevorzugt Daten der 
lokalen Umgebung eines aktuell betrachteten 
Iterationspunktes berucksichtigt werden. Ist eine zu 
geringe Menge lokalen Datenmaterials vorhanden, so wird 
dieses Datenmaterial durch gezielte Auswertung neuer 
Mefidaten (entsprechen Arbeitspunkten im Hinblick auf 
einen von den Zielf unktionen auf gespannten Parameterraum) 
geeignet erweitert (Experimental Design) . Die zu 
verwendenden Strategien zur gezielten Erweiterung des 
Datenmaterials hangen unter anderem von den zur Wahl 
stehenden lokalen Optimierungsverf ahren ab. Als 
Optimierungsverf ahren kommen in Frage: 

3.1. Simplexgestiitzte Optimierungsprozeduren (siehe 
[2]), die nur Zielf unktionswerte verwenden. 

3.2. Gradientengestutzte Optimierungsprozeduren. 

3.3. Prozeduren, die basierend auf der Idee der simplex 
gestutzten Algorithmen zusatzlich 
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Gradienteninf ormation verwenden. Mit den 
erweiterten Daten kann somit der nachste lokale 
Optimierungsschritt ausgefiihrt werden (Online- 
Op timie rung) . 



Durch Verwendung spezieller Kleinstquadratschat zer (etwa 
Kalmann-Filter, siehe [3]) mit AusreiAeridentif ikation konnen 
alle verwendeten Modelle online adaptiv verbessert werden, 

10 Spezielle Fragestellungen konnen zu 

Vektoroptimierungsproblemen fuhren. Das Ziel der 
Anlagensteuerung wird dabei durch eine vektorwertige Funktion 
F bestimmt: 



15 



F : R n x Q -+ R 



(40) , 




(41) 



mit einer zu wahlenden vektorwertigen Funktion g gemaft 



20 



g : R m x R n -> R d 



(42) , 



(r|, x) -> g(r|, x) 



(43*) . 



25 



Analog zu skalarwertigen Optimierungsproblemen konnen 
Vektoroptimierungsprobleme durch entsprechendes dreistufiges 
Vorgehen behandelt werden. 
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DEFINITION PROBLEMSPEZIFISCHER Z I EL FUNKT I ONEN 
(Fall d=l) 

In vorausgegangenen Abschnitt wurde die Struktur moglicher 
Zielf unktionen beschrieben. Fur konkrete Fragestellungen 
mussen konkrete Funktionen g definiert werden. Eine mogliche 
Wahl von g lautet wie folgt: 



wobei Gi und G2 positiv (semi-) definite Matrizen darstellen. 
Diese Funktion wird f olgendermafien interpretiert : 

1. Ein Vektor r\ beschreibt einen gewiinschten Anlagezustand. 

Gesucht ist ein Punkt x, fiir den die Anlage moglichst 
genau den gewiinschten Zustand annimmt. Dieser Vektor f| 

wird vom Anlagen-Operator (in Form von Abweichungen vom 
Ist-Zustand) vorgegeben, ist daneben aber auch 
automatisch vorgebbar. 

2. Es gibt Anlagenzustande, die unter alien Umstanden 
vermieden werden mussen. Daher ist man darauf angewiesen, 
Stellgrofcen x zu wahlen r die nicht zu weit von 
unproblematischen Stellgrolien entfernt sind. In der 
obigen Funktion wird eine spezielle unproblematische 
StellgrofJe x beriicksichtigt . 

3. Zusatzliche Restriktionen konnen durch Straf-Terme 
additiv in die Funktion g eingebracht werden. Ein 
Beispiel hierfur lautet wie folgt: 

Die Funktion g soli minimiert werden, wobei die 
Nebenbedingungen 



g : R m x R n -» R 



(44) , 



(ti, x) (r| - fO T Gi(Ti - i\) + (x - x) T G 2 (x - x) 



(45) , 



h(x) = 0 und k(x) < 0 
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zu berticksichtigen sind. Dies erfolgt insbesondere durch 
Minimierung des folgenden Terms: 

g + m • h 2 (x) + n 2 ' k 4 (x) . 

5 
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Pa tent anspruche 

1. Verfahren zum Entwurf eines technischen Systems, 

a) bei dem Mefidaten eines vorgegebenen Systems anhand 
eines Ersatzmodells beschrieben werden; 

b) bei dem eine Mafizahl fur die Gute des Ersatzmodells 
ermittelt wird; 

c) bei dem das Ersatzmodell dahingehend angepafit wird, 
dafi es eine moglichst hohe Mafizahl fiir die Gute 
auf weist ; 

d) bei dem die Mafizahl fur die Gute des Ersatzmodells 
dadurch verbessert wird, daft zusatzliche Mefidaten des 
vorgegebenen Systems ermittelt werden; 

e) bei dem das hinsichtlich seiner Gate angepafite 
Ersatzmodell zum Entwurf des technischen Systems 
eingesetzt wird. 

2. Verfahren nach Anspruch 1, 

bei dem das Ersatzmodell ein Regressionsmodell ist. 

3. Verfahren nach einem der vorhergehenden Anspruche, 

bei dem das Ersatzmodell fur eine weitere Erh5hung der 
Mafizahl der Gute verfeinert wird, indem mindestens ein 
vorgebbarer Teilbereich des Ersatzmodell lokal angepafit 
wird. 

4. Verfahren nach einem der vorhergehenden Anspruche, 

bei dem die zusatzlichen Mefidaten dort ermittelt werden, 
wo das Ersatzmodell nur unzureichend genau das 
vorgegebene System abbildet. 
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5. Verfahren nach einem der vorhergehenden Anspriiche, 

bei dem die Gute anhand einer quadratischen Abweichung 
der Mefidaten von den durch das Ersatzmodell bestimmten 
Daten ermittelt wird. 

5 

6. Verfahren nach einem der vorhergehenden Anspriiche, 



bei dem die Mefidaten nach deren Gute, bezogen auf deren 
Abweichung von den durch das Ersatzmodell bestimmten 
Daten, sortiert werden und eine vorgegeben Anzahl von n% 



7 . Verfahren nach Anspruch 6, 

bei dem die n% schlechtesten Mefidaten nicht aussortiert 
werden, wenn sie in einem zusammenhangenden Bereich 



8 . Verfahren nach einem der vorhergehenden Anspriiche, 
bei dem die Menge der Mefidaten im Rahmen einer 
Vorverarbeitung reduziert wird. 

20 

9. Verfahren nach Anspruch 8, 

bei dem die Vorverarbeitung eine Klasseneinteilung der 
Mefidaten umfafit. 

25 10. Verfahren nach einem der vorhergehenden Anspriiche, 
bei dem die mittels Entwurf gewonnenen Daten zur 
Steuerung einer technischen Anlage eingesetzt werden. 

11. Verfahren nach Anspruch 10, 
30 zur Online-Anpassung der Steuerung fur die technische 



10 



schlechtesten Mefidaten aussortiert werden. 



15 



liegen. 



Anlage . 
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12. Anordnung zum Entwurf eines technischen Systems , 

mit einer Prozessoreinheit , die derart eingerichtet ist, 
dafi 
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a) Melidaten eines vorgegebenen Systems anhand eines 
Ersatzmodells beschreibbar sind; 

b) eine Mafizahl fur die Gute des Ersatzmodells 
ermittelbar ist, indem die Mefidaten des vorgegebenen 
Systems mit durch das Ersatzmodell bestimmten Daten 
verglichen werden; 

c) aus der Maflzahl fur die Gute das Ersatzmodell 
dahingehend anpaiibar ist, dafi es eine moglichst hohe 
Gute aufweist; 

d) die Mafizahl fur die Giite des Ersatzmodells dadurch 
verbesserbar sind, dali zusatzliche Mefidaten des 
vorgegebenen Systems ermittelt werden; 

e) das hinsichtlich seiner Giite angepaftte Ersatzmodell 
zum Entwurf des technischen Systems einsetzbar ist. 
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